Asymptotic analysis of a Selberg-type integral via hypergeometrics 
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Abstract. We show how to determine the asymptotics of a certain Selberg-type integral by 
means of tools available in the theory of (generalised) hypergeometric series. This provides an 
alternative derivation of a result of Carre, Deneufchatel, Luque and Vivo [ar^iv: 1003.5996]. 



In [2], Carre, Deneufchatel, Luque and Vivo consider the Selberg-type integral 

and they determine its asymptotic behaviour when N,a,b all tend to infinity so that 
a ^ aiN and 6 ~ 6iA^, where ai and bi are given non- negative real numbers. The reader 
is referred to the introduction of [2] for information on motivation from random matrix 
theory connected to random scattering theory to investigate this question. 

It should be noted that So{a,b) is a Selberg integral, which can be evaluated in a 
product/quotient of gamma functions (cf. [3]). This being the case, the asymptotics of 
So{a,b) is easily determined by means of known asymptotic formulae for the Barnes G- 
function (see [7]). Thus, it suffices to consider the quotient 

J ^ Sk{a,b) 
S'o(a, b) 

(this quotient is denoted by Ik/N in [2]) and determine its asymptotic behaviour. By 
(now) classical identities in the theory of symmetric functions, it is shown in [2, Cor. II. 3] 
that 

1 V^/ ,,^fk-l\{N-t)kia + N 



-""'N-kl^^ H W ia + b + 2N-i-2)k ' 
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where {a)m is the usual Pochhammer symbol defined by {a)m = oc{a + 1) ■ ■ ■ (a + m — 1) 
for m > 1, and (q;)o = 1. In view of this rather simple expression, which, after substitution 
of a by aiN and of b by biN, is a rational function in A^, the problem of determining the 
asymptotic behaviour of the resulting expression looks innocent. However, inspection of 
the expression (1) shows that this is not so innocent at all: the summands are all of the 
order of magnitude A^'^ for large N, while (see Theorem and Corollary below), in fact, the 
expression in (1) has a finite limit as — )■ oo in the described situation. Hence, what 
happens here is a large amount of cancellation of leading asymptotic terms. This is what 
one has to get under control, and this constitutes the challenge in determining this limit. 

To solve the problem, Carre, Deneufchatel, Luque and Vivo develop a diff'erence operator 
calculus in [2], which is of interest in its own right. 

The purpose of the present note is to demonstrate that the tools for solving this problem 
are readily available in the literature on (generalised) hypergeometric series, and it is clear 
that they may serve as well to solve problems of similar type. 

Namely, using the standard hypergeometric notation 



pJ^q 



tti. 



dr 



E 



qjm 



where {ai)m, etc. are again Pochhammer symbols. Equation (1) becomes 



Jk 



(iV + l)fc-i(a + 7V-l)fc 
kl{2N + a + b-2)k 



-4-^3 



l-N,l-k,2-a-N,3-a-b-k-2N 
2-a-k- N,l-k-N,3-a-b-2N ' 



(2) 



The point of this conversion is that one may now look in the literature of hypergeomet- 
ric series, which provides a large arsenal of utilities, such as transformation formulae for 
hypergeometric series etc., and see whether one (or a combination of more) of these for- 
mulae enable us to convert (2) into a different form from which the asymptotics of Jf~ can 
be directly read off. As I show below, this is indeed the case. To the informed reader 
it will come as no surprise that the computations below have been carried out using the 
Mathematica package HYP [5] . 

The crucial observation to make is that the 4F3-series in (2) is not just an "arbitrary" 
hypergeometric series, but that it is a (—1) -balanced series (cf. [1, p. 140]), that is, that 
the sum of the lower parameters of the hypergeometric series minus the sum of the upper 
parameters is —1. Unfortunately, there are no transformation formulae available for ( — 1)- 
balanced series, but there are plenty of them available for 1-balanced series (or, as is 
common to say, for balanced series), such as, for example (see [4, Appendix (III. 16), Q 1 1]; 
this is T4302 in HYP), 



4-^3 



;i 



A, B, C, -n 
E,F,l + A + B + C-E-F-n 

(A), {E + F-A-B)n{E + F-A-C)r 



{E)r,{F)^{E + F-A-B-C)n 

-n,E-A,F-A,E + F-A-B-C 



X 4F3 



E + F-A-B,E + F-A-C,l-A-n' 



(3) 
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where n is a non-negative integer. That is, we are just by "2 off' in applying a transfor- 
mation formula. 

In order to relate series with (small) integral differences between the parameters, the 
hypergeometric literature offers contiguous relations. One of the simplest such relations is 
(this is C14 in HYP), 



4-^3 



A, B, C, D 
E,F,G 



BCD 
' EFG 



4^3 



A,S + 1,C + 1,L» + 1 
E+1,F+1,G'+1 



+ 4i^3 



A-1,S,C,L» 
E,F,G 



If we iterate this contiguous relation, then we arrive at 



4-^3 



A,S,C,L> 
E,F,G 



^ , {B)r{C)r{D)r 

{E)r{F)r{G)r^ ' 



+ 



r-1 

E 

s=0 



{B)AC)s{D) 



^"(i?j, (F), (G)/^' 



A,B + r,C + r,D + r 
E + r,F + r,G + r ' ^_ 

A-1,B + s,C + s,D + s 



E + s,F + s,G + s 



(4) 



We now apply (4) with r — k — 1 to the 4F3-series in (2). The reader should observe that 
thereby all the 4F3-series in the sum over s will be 0-balanced, while the single 4 Fa-series 
outside of the sum will reduce to 1 because one of the upper parameters is 1 — A; + r = 0. 
To be precise, after writing the result in compact form, application of (4) to (2) yields 



k-l 



(N + !),_! ja + N- 1), (1 - A:), (2 - « - AQ, (3 - g - b - A: - 2A0, 
^■^^ k\ {a + b + 2N - 2)k {2 - a - k - N)s {1 - k - N)s {3 - a - b - 2N), 



X 4F3 



-N, l-k + s,2-a-N + s,3-a-b-k-2N + s 
2-a-k-N + s,l-k-N + s,3-a-b-2N + s 



Continuing in the same spirit, if we apply the contiguous relation (4) a second time, this 
time with r = A; — s — 1, then this will produce a double sum in which all appearing 
4F3-series are balanced. To be precise, we obtain 



Jk 

k-l k-s-1 

= E E 



{N + l)fc_i (g + AT - l)fe (1 - k)s+t (2 - g - N)s+t {3-a-b-k- 2N)s+t 
s=o t'o k\{a + b + 2N -2)k{2- a- k-N)s+t{l-k- N)s+t{3-a-b-2N)s+t 



X 4F3 



b-k-2N + s + t,-l- N,2-a- N + s + t,l-k + s + t 
-b-2N + s + t,l-k-N + s + t,2-a-k-N + s + t 



Now the transformation formula (3) can be applied, and, after some simplification, the 
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resulting expression collapses to 

k-l k-s-1 



(a 



)k-s-t- 



.1 (1 - a - N)s+t+i (k- s- t)s+t is + t + 2)k-s-t-i 



s=0 t=0 



• 4-P3 



k\{2-a-b-2N)k 

l-k + s + t,k,a + b + N-2,a + N 
s + t + 2,a-l,a + b + 2N-l 



fc-i 

= E 

m=0 



(a - l)k-m-i (1 - a - N)m+i {k - m)m {m + l)k- 



4-^3 



1 



kl {2-a-b-2N)k 

-k + m + l,k,a + b + N -2,a + N 
m + 2,a-l,a + b + 2N-l 

Now the limit N,a,b 00, in which a ~ aiN and b ~ biN, can be safely performed. 
Hence, we have proven the following result. 

Theorem. The limit of the quantity as N,a,b ^ 00 such that a ~ aiN and b ~ biN 
is equal to 



\^ J\ai + bi+2j \ a, J 



m+l 



m=0 
k—m — 1 



k-m-l\ (A; + £-l)!(m + l)! / (ai + l)(ai + 61 + 1 



{k - l)!(m + e+iy. V ai(ai + &i + 2) 



(5) 



This theorem does not give the same expression as in [2, Theorem IV.4]. However, 
once one has found some expression, it is easy to pass to an equivalent one by using 
hypergeometric transformation and summation formulae. (Since, in the above derivation, 
there is plenty of room for variation — for example, the transformation formula (3) could 
be applied in different ways, or we could have applied a different transformation formula 
for balanced 4F3-scrics, sec [6, Sec. 4.3.5] — , it is conceivable that there is a direct route 
leading to the corollary below, using the above ideas. We did not pursue this, however; in 
particular, we did not do a systematic search of such an alternative route.) 

Corollary ([2, THEOREM IV.4]). The limit of the quantity Jk as N,a,b ^ 00 such 
that a ~ aiN and b ~ biN is equal to 

(ai + l)\ 



lim Jfc = iV(-l)^Y 

3=0 ^ 



{a, + b, + 2)'^+J ^ 



k 

i + j + 1 



Proof. In the expression on the right-hand side of (5) we interchange the summations 
over m and £, and we write the (now) inner sum over m in hypergeometric notation. This 
leads to 

k-i 



hm Jk = 5](-l) 



k-e-1 ^1 



.k-e-i 



(«i + 1)^'+-'- («i + hi + ly 



(A- + ^- 1)! 



£=0 



(ai + 61 + 2)k+i 



e\ {£ + iy.{k-e-iy. 

2,-k + e+l ai + l 



e + 2 



ai 
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To the 2-Pi-series we apply the transformation formula (see [6, (1.7.1.3), sum reversed at 
the right-hand side]; this is T2106 in HYP) 



A, -n 
C 



—n, 1 — C — n 1 — z 
l + A-C-n' ~ 



where n is a non-negative integer. Writing the resulting 2 i^i -series as a sum over i, this 
gives 



k-l k-l-1 

£=0 i=0 



(/c + £-l)! 



k-i-i-i fk £ 1 

i j{e-i)\e\{k-e-i)\{k-i-i){k-i) 



(ai + l)^-^ (ai + bi + l)^ 
(ai + 61 + 2)k+i 

Now we replace the term (ai + bi + 1)^ by Xlj=o (j) ("-'-)"' l*^! +^1 + 2)^"-', we make the sum 
over £ the inner-most sum, and finally write the sum over £ in hypergeometric notation. 
In this way, the last equation becomes 



k-lk-j-l . s 

hm = y y (-i)'=-^-M^ ^ — — 



(fc+j-1)! 



]\{k-3-\)\{k-%-\){k-%) 



(ai + 1) 



k—i 



(ai + 61 + 2) 



"2-^1 



-fc + z + j + l,A; + j ■ 



The 2-^1-series can be evaluated by the Chu-Vandermonde summation (see [6, (1.7.7); 
Appendix (III.4)]; this is S2101 in HYP) 



A, -n 

C ' 



1 



jC-A), 

{C)n 



where n is a non-negative integer. After finally replacing ihyk — i— j — 1, the claimed 
result follows immediately upon minor rearrangement of terms. □ 
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